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ABSTRACT 

We investigate the ellipticity of the point-spread function (PSF) produced by imaging an un- 
resolved source with a telescope, subject to the effects of atmospheric turbulence. It is important 
to quantify these effects in order to understand the errors in shape measurements of astronom- 
ical objects, such as those used to study weak gravitational lensing of field galaxies. The PSF 
modeling involves either a Fourier transform of the phase information in the pupil plane or a 
ray-tracing approach, which has the advantage of requiring fewer computations than the Fourier 
transform. Using a standard method, involving the Gaussian weighted second moments of in- 
tensity, we then calculate the ellipticity of the PSF patterns. We find significant ellipticity for 
the instantaneous patterns (up to more than 10%). Longer exposures, which we approximate by 
combining multiple (N) images from uncorrelated atmospheric realizations, yield progressively 
lower ellipticity (as l/y/N). We also verify that the measured ellipticity does not depend on the 
sampling interval in the pupil plane using the Fourier method. However, we find that the results 
using the ray-tracing technique do depend on the pupil sampling interval, representing a gradual 
breakdown of the geometric approximation at high spatial frequencies. Therefore, ray tracing 
is generally not an accurate method of modeling PSF ellipticity induced by atmospheric turbu- 
lence unless some additional procedure is implemented to correctly account for the effects of high 
spatial frequency aberrations. The Fourier method, however, can be used directly to accurately 
model PSF ellipticity, which can give insights into errors in the statistics of field galaxy shapes 
used in studies of weak gravitational lensing. 

Subject headings: atmospheric effects - gravitational lensing 



Introduction 



Statistical analyses of w eak gravitational lens 



ing of field galaxies (e.g. , Wittman et al.1 2000: 



Van Waerbeke et al.|[200q iBacon et alB oOO) 



are 



being used as probes of cosmology and are ex- 
pected to provide some of the strongest cosmo- 
logical tests in future, large astronomical survey 
projects, such as the Large Synoptic Surv e y Tel e- 



scope (LSST; e.g., iTvson fc Angell 120011 : iTvson 



1 University of California, Department of Physics, 1 
Shields Ave, Davis, CA 95616 

institute for Geophysics and Planetary Physics, LLNL, 
L-413, 7000 East Avenue, Livermore, CA 94550 

3 Lawrence Livermore National Laboratory, 7000 East 
Avenue, Livermore, CA 94550 

4 University of Washington, Department of Physics, Box 
35160, Seattle, WA 98195 



2002). These surveys will allow for the precise 
determination of various cosmological parameters, 
such as the matter density distribution 12 m , the 
cosmological constant 57a, the equation of state w 
of the dark energy, and its time derivative. This 
is done by accurately analyzing large numbers of 
background galaxies as their shapes are sheared 
by intervening large-scale structure through weak 
gravitational lensing, the results of which are then 
combined with, e.g., the very accurate measure- 
ments of the cosmic microwave background radia- 
tion by the WMAP satellite. 

A critical part in these analyses is the accu- 
racy to which one can measure and correct the 
shape of the Point Spread Function (PSF) as 
it va ries across the detector (see, e.g., Irloekstral 



2004). This PSF anisotropy is largely induced by 
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the atmosphere (e.g., IWittmanl |2005| ). and can- 
not easily be modeled without incorporating an 
explicit atmosphere. It is possible to mimic the 
effects of the atmosphere by convolving either ar- 
tificially generated, or high resolution HST im- 
ages with a suitable PSF (see, e.g., Heymans et 
al. 2006a for the former, and Bacon et al. 2001 
for the latter approach), but that still does not in- 
clude effects of PSF anisotropy. We therefore set 
out to model the behavior of the PSF as it gets 
folded through a realistic atmosphere and tele- 
scope system. Since ray-tracing methods are com- 
monly used to simulate the shearing sign al of weak 



lensing by large scale structure (e.g., Jain et aD 



tensing 

120001: IVale fc White! 120031 : iHevmans et al.ll2006ri) . 
we include both ray-tracing and Fourier trans- 
form methods to calculate what the PSF should 
be based on the phase and intensity information 
in the pupil plane of the telescope. The ray- 
tracing method has the advantage of computa- 
tional speed compared with the Fourier transform 
method. However, the Fourier transform method 
correctly treats the effects of interference, which 
are ignored by the geometric ray-tracing method. 
Because of its relative computational efficiency, it 
is useful to understand whether the ray tracing 
approach gives adequate PSF ellipticity informa- 
tion. More importantly, it is crucial to understand 
the elliptical properties of the PSF in order to en- 
able quantitative analysis of the errors in statisti- 
cal studies of field galaxy shapes. 

Our modeled observational setup is character- 
ized as follows: a generic 8 m-class telescope, a 
turbulent atmosphere with a Kolmogorov power 
spectrum, and a single, on-axis point-source lo- 
cated at infinity. The ellipticity and its direction 
are assumed to be representative for a single re- 
gion over which these quantities do not vary. The 
angular extent of these regions can be, depending 
on o bserving conditions, l arger than an arcminutc 
(e.g. JAsztalos et al.ll2006l ). Since the typical LSST 
exposure will be 15 seconds, we also investigate the 
time dependence of the atmospherically induced 
ellipticities, and how it imposes limits on the abil- 
ity to measure them. 

1.1. Layout of the paper 

The paper is organized as follows. In § [2] and 
§ [U we discuss the atmospheric simulations and 
the methods applied to calculate the PSF elliptic- 



ity. Then in § [4] we describe the results for both 
the ray-tracing method (RTM) and the Fourier 
transform method (FTM). The main difference be- 
tween these two is the way one models image for- 
mation by the telescope. The RTM traces the ge- 
ometric path of rays from the pupil plane onto the 
focal plane, whereas the FTM applies a Fourier 
transform to the complex electromagnetic field in 
the pupil plane in order to calculate the resulting 
image. The latter correctly incorporates effects 
of interference, unlike the RTM. The geometric 
approximation of the RTM has implications for 
shape measurements in the image plane, which we 
quantify in this paper. 

The results are subdivided into the effects of 
pupil plane sampling (§ 14.1)1 . the variation of el- 
lipticity as a function of exposure time and the 
presence of wind (§ 14. 3|) . and seeing (S I4.4p . The 
latter is approximated by using varying ratios of D 
/ r (where D is the aperture size of the telescope, 
and r Q is the coherence length of the atmosphere) . 
Longer exposure times are simulated by increas- 
ing the number of completely independent phase 
screens (which is a function of aperture diameter 
and wind-speed). In order to check the accuracy of 
this approximation, we also evaluate a model that 
includes intermediate phase screens, i.e., screens 
that are not completely decorrelated from the pre- 
vious one, but are translated by a small fraction of 
the aperture size along the wind direction. In each 
of these sections we investigate the differences be- 
tween the RTM and FTM, which are summarized 
in§El 

2. Simulations and Methods 

We generated random phases for the electro- 
magnetic field in the pupil of the telescope using 
Kolmogorov statistics to represent the effects of at- 
mospheric turbulence. The inner turbulence scale 
of the simulations is set by the pixel size used in 
the simulations, to a fraction of r Q . The outer 
scale of the turbulence has been fixed to a value 
larger than the simulation box (> 800 m for the 
large screens in § 14.3) for instance), so that ef- 
fectively the simulations see a Kolmogorov turbu- 
lence spectrum with an infinite outer scale. As- 
suming a constant magnitude in the circular pupil, 
we propagated the field to the focal plane using a 
Fourier transform, and by squaring the resulting 
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Fig. 1. — Representative point-spread functions, 
using the Fourier method (left two images), and 
the ray-tracing method (right two images). The 
top images are for an instantaneous realization of 
the atmosphere, whereas the bottom images are 
the means for 256 such instances. Notice that 
the top left Fourier image displays a prominent 
speckle pattern due to interference. This pattern 
gets washed out over time. 

magnitude, created a representation of the focal 
plane image intensity. These images contain dis- 
tinct speckle patterns for different realizations of 
the atmosphere. 

The ray tracing algorithm used an idealized 
model of a reflecting telescope. A Kolmogorov 
phase screen was placed in the aperture of the 
telescope and a uniform distribution of rays in the 
pupil was assumed. The x and y derivatives of 
the phase screen were used to determine the atmo- 
spherically induced deviations of the rays as they 
propagated toward the focus of the telescope. At 
the focus, the incoming rays are mapped onto a 
fiducial 2D detector grid which determines the in- 
tensity distribution in the focal plane. 

3. Calculating Ellipticity 

We calculate the ellipticity of an object in 
the pupil plane as follows. Assume we have an 
image for which the pixel coordinates are given 
by x = l..N,y = 1..N; the intensity in each 
pixel is given by I(x,y); the seeing is given by a 
(FWHM=2. 355cr for a Gaussian), and the central 



point source is located at (x c , y c ). We then define 
Gaussian weights as follows: 



w[x,y) — — p c > . — — e 



V(W) 



(1) 



and the Gaussian weighted moments: 



N N 

x,y — l x,y — l 

N N 

s x = ^2 ( xI ( x ,y) w ( x ,y)), Syy = (y 2l ( x ,y) w ( x ,y)) 

x,y — l x,y — l 

N N 

Sy = ^2 (yl(x,y)w(x,y)), S xy = ^ (xyl(x,y)w(x,y)) 

x,y — l x,y — \ 



This allows us to define the following quantities 
(all weighted by I(x,y)w(x,y)): 



(2) 



x c — S x / S w , y c — Sy/ S w 




SyyS W SySy (3) 




From these, one can calculate the image Gaussian 
scale length a and the ellipticity e: 



(7 = ^/rxx + ryy 

1/2 

({rxx — ryy) 2 + (2rxy) 2 \ 



(rxx + ryy) 2 
rxx — ryy Irxy 
rxx + ryy ' 



c-2 



rxx + ryy 



(4) 



This i s the method used by Kaise r. Squires, fc Broadhurst 
(|l995l KSB) - see also iHevmans et al.l (|2006al ) for 
an overview of the different weak lensing pipelines. 
The weights in eqn. 1 are dependent on initial val- 
ues of (J, and the source centroid (x c ,y c ). How- 
ever, one can iterate from initial guesses for the 
unknowns. This algorithm converges quickly (typ- 
ically within a few steps), and yields values for a 
and the source position that are readily verifiable. 
We terminate the iterations when the changes in 
a are less than 0.01 pixel. 
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Fig. 2. — Ellipticity as function of the number of independent phase screens and sampling rates for D / 
r = 40 (r Q = 21cm). The left panel shows the results for ray-tracing, with the pupil-plane sampling rates 
color-coded as: r Q / 2 green, r Q / 4 blue, r D / 5 cyan, r D / 8 purple, and r D / 16 yellow. The panel on the 
right shows the same results (except r Q / 5) for the Fourier method. The latter method clearly illustrates the 
expected behavior: ellipticity should be independent of pupil plane sampling rates beyond Nyquist rates. 
The progressive lowering of the curves for higher samping rates in the left panel, therefore, is unphysical. 
The dashed line in the right panel shows a 1/VN decline. 



4. Results 

Most of the figures in this paper show the ellip- 
ticity behavior as a function of the number N of 
independent phase screens. A single phase screen 
(N = 1) therefore represents the instantaneous el- 
lipticity of a particular representation of the at- 
mosphere (see Fig. [I]). In all of our subsequent 
discussions we assume that the telescope is per- 
fect, i.e., it does not induce image aberrations. 

Whether we apply the RTM or the FTM ap- 
proach, we first create stacks of 500 completely 
uncorrelated instances of the atmosphere (actu- 
ally phase-screens in the pupil plane). The result- 
ing focal plane images are then either ray-traced 
or calculated using an appropriate Fourier trans- 
form. After this, we randomly select N frames 
out of the 500, which are then stacked, averaged, 
and have their ellipticity calculated. We repeat 
this 1024/A times. Since there are not 1024/AT 
fully independent stacks present (for N > 2), 
some smoothing occurs, especially for the larger 
N stacks. The figures show the mean ellipticity 
for these 1024/iV stacks, and the error-bars on the 
means are approximated by the rms of the distri- 
bution divided by the square root of the number 



of stacks (1024/iV). 

4.1. Constant ratio D / r Q , varying sam- 
pling in pupil plane 

The value of D / r Q has been fixed at 40 (r Q = 
21 cm) for an assumed 8.4 m pupil diameter with 
a central obscuratior0 of 5.4 m, while the (phase) 
sampling rate in the pupil plane increases from r / 
2 to r Q / 16. The ellipticity e as function of the 
number of independent phase screens N is shown 
in Fig. [2] The left panel illustrates the results for 
the RTM method. While each individual sampling 
rate falls off as 1/yN, they are offset in ellipticity 
as the sampling rate increases. An increase from 
r / 2 to r Q / 16 more than halves the measured 
ellipticity for a given number of independent phase 
screens. This is clearly not physical. 

The following example may illuminate this 
behavior. Assume, for instance, a simple one- 
dimensional cosine phase wave with frequency k. 
The Fourier transform of this function produces 
two delta functions located at ±fc on the x-axis. 
However, since ray-tracing uses the derivatives 
of the phase to calculate where the rays will go, 



1 These are the current parameters for the LSST design 



it produces a spread of points due to the range 
in derivatives. An increase in the sampling rate 
will increase the likelihood high angle rays will 
be modeled. Furthermore, the derivative range 
becomes larger for higher spatial frequencies (it 
goes from — k to +fc) resulting in progressively 
more aberrant rays, regardless of the sampling 
rate. While perhaps an extreme example, it does 
underline the fact that 2D ray-tracing will pro- 
duce a more spread-out image due to rays being 
deflected into unphysical angles. This broadening 
of the image then results in a lowered ellipticity as 
it decreases the local asymmetry (remember that 
the ellipticity contribution of a point is weighted 
by its distance, see Eqn. 1). 

The FTM, on the other hand, does display the 
correct behavior (for the exact same sets of phase 
screens): the ellipticity is independent of the sam- 
pling rate (provided it is at least r Q / 2). Except 
for the sampling dependency, both methods ex- 
hibit the following characteristics. First, elliptic- 
ities decrease linearly (in log-log) as the number 
of frames increases. The slope is consistent with 
a 1/yN decline (a = —0.5), as indicated by the 
dashed line in the right panel. And second, ellip- 
ticities of individual speckle images (N = 1) are 
-9% for D / r c = 40 (r Q = 21 cm). 

It should also be noted that there are ~ 20 — 40 
independent instances of the atmosphere for an 
8.4 meter telescope aperture with wind-speeds of 
~ 10 — 20m/s (typical, turbulence- weighted values 
for many astronomical sites) and a 15 second expo- 
sure. Therefore, these simulations predict that the 
raw ellipticity of a 15-second exposure for a point 
source image through an 8.4 meter telescope with 
D / r Q = 40 is ~ 2%. We explore the effects of 
wind in more detail in § 14.31 

4.2. Binning in the image plane 

Obviously, no astronomical instrument de- 
signed for seeing-limited observations will sample 
the PSF at the Nyquist interval for the telescope 
diffraction pattern (we measure a ~ 60 pixels), 
so the next step is to see what happens to the 
ellipticitics if one progressively robins the images 
of Fig. [2 The results are listed in Table 1. The 
ellipticities are for a stack of 32 random frames 
with D / r Q = 40. This stack gets increasingly 
rebinned down to scales where the PSF is barely 
resolved (a ~ 1). 
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Fig. 3. — Effect of wind on ellipticity behavior for 
D / r = 40 (r = 21 cm). The blue dashed curve 
is for phase screens which are translated across 
the aperture, and the red curve is for independent 
instances of the atmosphere. Both these curves are 
calculated with the Fourier method. The wind- 
speed v is needed to convert the number of frames 
N into an elapsed time t (— ND/v, with D — 8.4 
m, and v in units of m s _1 ). Therefore, a typical 
15 s LSST exposure, with a wind-speed of 10 m 
s _1 , contains 18 frames. 

Given that the values of e do not change sig- 
nificantly over a large range of binning, it is clear 
that the ellipticity measurements are robust and 
do not depend on the pixel scale. For compari- 
son, the rms spread in the value of e for distinct 
random stacks of 32 images is about 60% of the 
mean value of e, whereas the listed relative range 
in column 3 is only ~ 5% under rebinning. 

4.3. Atmospheric model with wind 

So far we have only considered the behavior of 
ellipticity as function of the number of uncorre- 
cted instances of the atmosphere. As described 
in § 14.11 the evolution of the PSF with increas- 
ing exposure time can be estimated using the re- 
sults from these uncorrelated screens by associat- 
ing each screen with a unit of time equal to the 
aperture diameter divided by the wind speed, i.e., 
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Table 1 

Pixel binning effects on ellipticity 



Binning n 


a n [pix] 




Ratio e„/ei [%] 


1 


61.82 


0.02505 




2 


30.92 


0.02507 


100.1 


4 


15.47 


0.02504 


100.0 


8 


7.74 


0.02513 


100.3 


1G 


3.89 


0.02492 


99.5 


32 


1.97 


0.02393 


95.5 


64 


1.04 


0.02443 


97.5 



Note. — The ellipticity is calculated on a random stack 
of 32 speckle patterns, with a sampling of r Q / 2, and a 
ratio of D / r Q = 40 (r = 21 cm). There is no significant 
dependency on pixel size. Note that a ~ 1 to 2 are typical 
astronomical seeing disk sampling ratios. 



the time it would take for the wind to translate a 
screen completely out of the aperture. However, 
a more realistic treatment of the PSF evolution 
involves a more continuous translation of a Kol- 
mogorov phase screen across the telescope aper- 
ture. 

For this purpose, we generated 3 large phase 
screens which contain 95 aperture-clearings each. 
The translation offset0 is such that an aperture- 
clearing takes 20 steps; from each phase screen 
we therefore generate 1900 pupil images using the 
FTM. The ellipticity is then calculated on com- 
bined stacks of 20 pupil images, yielding one value 
per aperture-clearing. The results are plotted in 
Fig. El The blue dashed line shows the mean ellip- 
ticity values for all 3 x 95/iV independent aperture- 
clearings ("frames"), with as error-bars the error 
in this mean. At N — 5, for instance, we cal- 
culated the mean of all uncorrelated instances of 
5 x 20 consecutive pupil images. On the other 
hand, the red curve shows the ellipticity behav- 
ior of the sum of N uncorrelated instances of the 



2 This offset should not be confused with the wind-speed v. 
All we want to make sure is that we have enough numerical 
resolution (hence the 20 steps) as the atmosphere translates 
across the aperture. How long it takes for the atmosphere 
to clear an aperture does not matter for this calculation. 
See also the caption to Fig. [3] 



atmosphere, using the same pupil plane sampling 
and value of D / r Q = 40 (r Q = 21 cm). 

A few things stand out. Below about 10 frames 
or so, the curves behave differently. In case of 
the red line, the ellipticity jump from N = 1 to 
N = 2 is due to the fact that there is no "image 
motion" in the N = 1 case, whereas for N = 2, 
two pupil images have been combined with differ- 
ent PSF centroids. This raises the ellipticity be- 
yond what is there in a single PSF (one would need 
to shift-and-add to remove this effect). It subse- 
quently takes a few more co-added frames for this 
centroid-offset effect to cancel out (on average the 
PSF has to align with the optical axis since we put 
the source there). 

The blue line does not suffer from this centroid- 
offset problem since the phase-screens are contin- 
uous (only l/20 th gets shifted out between pupil 
images) and the PSF centroid cannot move around 
discontinuous as a consequence. However, we do 
see another effect present in the blue curve. Be- 
cause the Kolmogorov phase screen will generally 
have low-spatial-frequency correlations that are 
larger than the telescope aperture, the ellipticity of 
the PSF is expected to decrease more slowly with 
increasing N than for the discontinuous model 
using multiple independent phase screens. After 
some time (or equivalently, for larger values of TV), 
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Fig. 4. — Speckle patterns for individual phase 
screens at various ratios of D / r Q (r Q = 21, 42, and 
84 cm, respectively). Note that with increasing 
values of r (left to right), the number of speckles 
decreases, while their intensities go up. Also, the 
size of the pattern decreases with increasing r . 

these low-spatial-frequency correlations disappear 
and the slope of the blue curve steepens to that 
of the red curve. This appears to be happening 
between N = 10 and N = 20. 

The main result of this exercise, however, is the 
confirmation that our method of using uncorre- 
cted instances timed at a rate equal to an aper- 
ture clearing time is a valid approximation for the 
ellipticity behavior in a statistical sense, as long 
as we are in the long-exposure, N > 20 domain. 

4.4. Varying r Q 

We also investigated the effect of increasing the 
value of r Q . The expectation is that for larger val- 
ues of r D (i.e., better seeing conditions) the num- 
ber of speckles goes down (no atmosphere = no 
speckles, diffraction pattern only) while their in- 
dividual brightness goes up (due to the conserva- 
tion of flux) . This is easily verified in the individ- 
ual speckle patterns (see Fig. [4|. It is not clear, 
however, what the behavior of the RTM method 
with respect to the FTM will be. In § 14.11 we no- 
ticed that the RTM method significantly under- 
estimates the actual ellipticity depending on the 
sampling rate. If the relative offsets are constant 
then one might be able to come up with a par- 
ticular sampling rate for the ray-tracing case that 
best matches the actual ellipticities (for the D / 
r Q = 40 case the best matching sampling looks to 
be about r Q / 3, see Fig. [5]). 

The results are presented in Fig. [5l with the 
solid lines representing the FTM results, and the 
dashed lines are for the RTM method. If we fo- 
cus on the solid lines first, it is clear that the el- 
lipticities increase with increasing values of r Q at 




1 10 

Frames (N) 

Fig. 5. — Measured ellipticities e (see eqn. 4) for 
different values of r , with values from top to bot- 
tom of 42 (red, D / r o =20), 21 (green, D / r Q = 
40), and 10.5 cm (blue, D / r = 80), respectively. 
The solid line curves have been calculated using 
the Fourier method, and have a fixed pupil plane 
sampling rate of r Q / 4. The dashed lines are cal- 
culated using ray-tracing, and are color-coded and 
sampled similarly. The ellipticity for a given num- 
ber of independent phase screens depends on the 
size of r Q : large values of r Q have larger ellipticities. 
Also note that low number statistics are affecting 
the data-points toward large frame counts causing 
the curves to cross eachother. 

the same number of stacked frames. This can be 
qualitatively understood in terms of the decreas- 
ing number of speckles distributing themselves in 
a progressively less circular pattern due to the 
smaller number statistics. 

Based on this plot, it is also apparent that, 
even though the RTM underestimates the elliptic- 
ity compared to the FTM, it does so more or less 
independently of the value of r Q (which is a proxy 
for seeing). This might open up the possibility 
that one either selects a computationally efficient 
RTM sampling rate (say, r Q / 2) and apply an ap- 
propriate (fixed) correction factor to the elliptic- 
ity results, or adjust the sampling rate such that 
the RTM and FTM results are in good agreement 
(~r / 3). Another approach that is under inves- 
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tigation (G. Jernigan, private communication) is 
to roll off the atmospheric power spectra at high 
spatial frequencies. Further study is needed to as- 
sess the accuracy of any of these approaches. For 
instance, the numerical correction factors derived 
from the results shown in Fig. 5, for the particular 
sampling ratio of r Q / 4, are 1.32, 1.15, and 1.21, 
for r Q = 42 cm, r Q = 21 cm, and r Q = 10.5 cm, re- 
spectively. Whether this variation is due solely to 
the statistical errors in our modeling, or includes a 
systematic dependence on r Q is not known. Clearly 
if one requires accurate modeled ellipticities, then 
the computationally more expensive FTM is cur- 
rently preferred. 

5. Conclusions 

Based on this analysis, we reach the following 
conclusions: 

1. Instantaneous speckle patterns have elliptic- 
ities of -10% for D / r Q = 40 (r Q = 21 cm). 

2. Co-adding multiple patterns results in a 
linearly decreasing ellipticity (on a log- 
log plot), consistent with a s/N slope of 
a = -0.5. 

3. Modeling phase screen transport across the 
aperture (i.e., wind) does not change the el- 
lipticity results obtained from adding uncor- 
rected phase screens in the limit of long ex- 
posures (N > 20 aperture clearings). 

4. Ellipticity values arc robust over a large 
range of pixel binning. 

5. We expect no ellipticity dependency on sam- 
pling density in the pupil plane (above sam- 
pling of r Q / 2). This is confirmed for the 
Fourier method, but not for ray-tracing. The 
latter has a strong dependency on sampling 
rate, in the sense that the higher the sam- 
pling rate, the lower the resulting ellipticity. 
This can be understood as the result of a 
breakdown in the geometric approximation 
for high spatial frequency aberrations. 

6. Ellipticities grow (for a given N) as the value 
of r Q increases. However, since the average 
size of the PSF goes down as r Q increases 
(see Fig. 2]), the ability to measure precise 



ellipticities actually improves (for a given re- 
solved object). 

In summary, the effects of interference must be in- 
cluded in order to comprehensively model point- 
source ellipticities induced by the atmosphere. 
Therefore, care has to be taken that the geomet- 
ric optics approximation to image formation by 
the telescope (i.e., ray tracing) produces the same 
modeling results, as this is not true in general. 

The authors like to thank Don Phillion for gen- 
erating the large phase screens of § 14. 3i and the 
anonymous referee for helpful comments. This 
work was performed under the auspices of the 
U.S. Department of Energy, National Nuclear Se- 
curity Administration by the University of Cali- 
fornia, Lawrence Livermore National Laboratory 
under contract No. W-7405-Eng-48. 
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